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Many inorganic pigments contain heavy met- 
als hazardous to health and environment. Much 
attention has been devoted to the quest for non- 
toxic alternatives based on rare-earth elements. 
The computation of colors from first principles is 
a challenge to electronic structure methods how- 
ever, especially for materials with localized /- 
orbitals. Here, starting from atomic positions 
only, we compute the color of the red pigment 
cerium fiuorosulfide CeSF, as well as of mercury 
sulfide HgS (classic 'vermilion'). Our methodol- 
ogy employs many-body theories to compute the 
optical absorption, combined with an intermedi- 
ate length-scale modelization to assess how col- 
oration depends on film thickness, pigment con- 
centration and granularity. We introduce a quan- 
titative criterion for the performance of a pig- 
ment. While for HgS this criterion is satisfied due 
to large transition matrix elements between wide 
bands, CeSF presents an alternative paradigm: 
the bright red color is shown to stem from the 
combined effect of the quasi two-dimensionality 
and the localized nature of 4/ states. Our work 
demonstrates the power of modern computational 
methods, with implications for the theoretical de- 
sign of materials with specific optical properties. 

Light propagating inside a heterogeneous solid experi- 
ences (i) absorption, and (ii) scattering. The light that is 
not absorbed is diffusely reflected and is responsible for 
the perceived color. The visual appearance of a material 
is hence determined by selective absorption of light and 
sufficient (back) scattering. For a material to be e.g. a lu- 
minous red pigment, two criteria must thus be satisfied. 
First, its absorption edge should be located at the appro- 
priate energy (~ 2.1 eV) so that the red component of the 
visible spectrum is not absorbed. Second, the absorption 
edge should be sharp, so that most other photons within 
the visible range (green, blue) are absorbed. 

The computation of these effects from first principles 
is faced with three fundamental difficulties. First, in 
view of the sensitivity of the human eye, the optical gap 
must be obtained with a precision of at least 100 meV. 
Conventional electronic structure methods yield a well- 



documented underestimation of the gap of conventional 
semiconductors. Second, the localized 4f states, which 
play a crucial role in optical properties of rare-earth based 
pigments [H-Q, arepoorly described by standard density- 
functional theory [j] or even G W approaches (loj . Third, 
a realistic assessment of the coloration of a pigment must 
take into account scattering properties depending on con- 
centration, granularity and film thickness. Ab initio sim- 
ulations so far have not ventured beyond calculating the 
optical conductivity of infinite bulk samples (see however 
Ref. HQ for organic molecules). In this article we ad- 
dress all these issues and develop a general methodology 
for the prediction of the color of narrow-band materials. 

We investigate cerium fiuorosulfide (CeSF), a typical 
example of the new class of rare-earth pigments 0, [ii} . It 
crystallizes in the layered ThCr2Si2 structure sketched in 
Fig. Q] (which is incidentally also that of the recently dis- 
covered iron-based superconductors [ltj). The same fig- 
ure displays the momentum-dependent many-body spec- 
tral function Ak(u)) that encodes the excitation energies 
associated with the addition or removal of an electron 
into the many-body ground-state (see Supporting Infor- 
mation). The localized Ce-4/ states form quasi-atomic 
multiplcts that hybridise weakly with the rest of the solid. 
From the panels (c)-(e), displaying the specific orbital 
character of each electronic state, it is apparent that the 
top of the valence band has dominantly S-3p character. 
The occupied Ce-4/ states are located at higher binding 
energies, near the center of the S-3p bands. The lowest 
unoccupied states are found to be the almost localized 
empty Ce-4/ states, with the dispersing Ce-5d conduc- 
tion band lying higher in energy. 

The calculated optical gap A = 2.14 cV (578 nm) is 
consistent with the red color of this material. Further- 
more, our calculation reveals that the absorption edge is 
associated with the optical transition between S-3p states 
and Ce-4/ states. This is at variance with the standard 
lore associating the coloration of cerium-based pigments 
to the intra-atomic Ce-4/ to Ce-5d transition. Indeed, 
such a conventional assignment was proposed for CeSF in 
Ref. 0, and X-ray photoemission spectra (XPS) were in- 
terpreted along the same lines in Ref. [H|]. m figure [3 we 
show that our results are nonetheless fully consistent with 
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FIG. 1: Cerium fluorosulfide CeSF. a) Crystal structure. The red, yellow and grey spheres represent Ce, S, and F, 
respectively, b). The total many-body momentum-dependent spectral function Ak(w) along high symmetry lines, c-e). The 
corresponding partial spectral functions for Ce 4/ (c), S 3p(d), and Ce 5rf(e). f). The total momentum- integrated spectral 
function (black curve) as well as the partial momentum-integrated spectral functions for Ce 4/ (red), S 3p (blue) and Ce 5d 
(green) . 
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FIG. 2: X-ray photoemission spectroscopy (XPS). The 

many-body spectral momentum-integrated function (black 
curve) of CeSF compared to experimental XPS spectra [lit] 
(thick green curve). In order to simulate experimental reso- 
lution the theoretical spectral function has been convoluted 
with a Gaussian of full width at half maximum of 0.3 eV. 
As a comparison we show (red curve) the XPS spectrum of 
the alternative intra-atomic scenario previously proposed [See 
the text and Supporting Information for a discussion]. The 
shaded regions of the corresponding colors are the contribu- 
tion from Ce 4f states in each case. 



assignment of the observed spectral peaks: according to 
our calculations, the peak at -3 eV is due to contributions 
from localized Ce 4/ states while the shoulder at -1.5eV 
is mainly formed by S-3p states that are pushed upwards 
in energy by the hybridization with Ce-4/ states. The 
previously proposed 'intra-atomic' scenario can also be 
simulated theoretically (see Supporting Information) and 
is shown in Fig. [2] to be inconsistent with the measured 
spectra. 

In contrast to CeSF, mercury sulfide cv-HgS (also 
known as cinnabar red or "vermilion") has been used 
as a red pigment since antiquity [l2[ . It is a conventional 
semiconductor with a simple hexagonal structure. Al- 
though the additional complication of localized states is 
not present here, an accurate determination of the gap 
requires the use of G W calculations (see Supporting In- 
formation for details). The onset of absorption is due to 
transitions between broad and strongly hybridizing bands 
of mainly S-3p and Hg-6s character. Given the qualita- 
tively different nature of the optical transitions involved, 
HgS and CeSF are good pigments for entirely different 
reasons, as explained in detail below. 

The absorption properties of these two compounds can 
be derived from the frequency-dependent optical conduc- 
tivity, Re<r(a;). Using linear response theory, and ne- 
glecting excitonic effects (see the Supporting Information 
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FIG. 3: Optical conductivity and absorption coefficient of CeSF and HgS. (A) Polarization averaged optical conduc- 
tivity. The inset is a zoom into the energy region of the onset of absorption. (B) Macroscopic absorption coefficient K. The 
horizontal dashed lines are the quality criterion ([3| as defined in the text computed for f3 — 50/mm. 



for a discussion), it can be expressed as (see e.g. [13. 14|) 

Re d a (n) = |^ E j d JM-f(^) 
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x tr|Ak(w + 0.)hvk, a Ak(uj)hv^ a | (1) 

where V is the unit-cell volume, a labels the polarization, 
Vk. a are the optical transition matrix elements and the 
Fermi functions f(w) restrict transitions to take place be- 
tween occupied and empty states. The absorption is de- 
scribed by the macroscopic absorption coefficient K(u>): 



K(u) = —k(cj) = —ImJ 1 + —a(u) (2) 

C e C e V W 

where k(u>) is the imaginary part of the refractive index, 
and c e the speed of light. 

The calculated optical conductivity and absorption 
coefficient of CeSF and a-HgS arc displayed in Fig. [3] 
(a),(b). While marked differences exist between the two 
compounds on a broad energy range, the magnitude of 
a and K near the onset of absorption (see the inset 
of figure [3ta) and figure EJb)) is actually similar, with 
o~ — 50 O^cm" 1 and K ~ 20nmr 1 for ui = A + 0.1 eV. 
The physical origin of the spectral weight just above the 
absorption edge (which in turn determines K) is however 
fundamentally different for HgS and CeSF. Indeed we will 
now show that the large absorption in HgS is due to the 
strength of optical transition matrix elements, while for 
CeSF it is due to the large density of localized Ce-4/ as 
well as weakly dispersing S-3p states. 

To substantiate this claim, we note that color is de- 
termined by a fairly small frequency range above the 
absorption edge (see also below), and, hence, we focus 
on A < hbj < A + hSu c , with HSuj c w 0.2 eV. In or- 
der to disentangle density of states effects from tran- 
sition probabilities, we compute the ratio between the 
optical conductivity, equation ([T}, and the joint den- 
sity of states N(u)) [3l[ at the energy uj = A + h5u c . 



We find that a/N = 2.5 (h/(a m e )) 2 for HgS, while 
a/N = 0.86 (h/(aom e )) 2 is almost three times smaller 
for CeSF (ao is the Bohr radius and m e the mass of the 
electron). This proves that transition matrix elements 
dominate for HgS, while density of states effects domi- 
nate for CeSF. The reason for this is the weak dispersion 
of the Ce-4/ states (expected from their localized char- 
acter), but also of the S-3p states at the top of the va- 
lence band in CeSF. Interestingly, the weak dispersion of 
the S-3p states along the c-axis is due to the quasi two- 
dimensional environment of the sulfur atoms, which are 
located within layers parallel to those containing cerium 
atoms (Fig. QJa)). As detailed in the Supporting Infor- 
mation, the joint density of states N(uj) and the absorp- 
tion K(uS) display a discontinuous increase at threshold 
N oc 9(fud — A) for a strictly two-dimensional disper- 
sion. In contrast, a three-dimensional parabolic disper- 
sion (mimicking e.g. the 3p or 6s bands in HgS) yields 
a less sharp frequency-dependence oc (tvjj — A) 1 / 2 above 
the absorption edge. 

The macroscopic quantities a and K describe the ab- 
sorption properties of a perfectly crystalline bulk solid. 
Pigments are used however (in e.g. paints) in the form of 
small particles embedded into a transparent matrix. In 
this context, light propagation depends on the morphol- 
ogy of the sample, and a multiple scattering problem has 
to be solved to determine the diffuse reflectance RjuS] (in 
contrast to the simpler specular reflectivity [29|, H(|). A 
commonly used [35J approach is the effective-medium de- 
scription due to Kubelka and Munk (KM) [HHi]. Thc 
KM model treats the propagation of light through a ho- 
mogeneous layer of matter with pigment concentration 
c, that absorbs light with rate cK(lo) and backscatters 
with rate /3(ui) per unit length. The quantity /? effec- 
tively contains all the information on the microscopic 
structure of the sample. As shown in the Supporting 
Information, the energy dependence of /3 can be ne- 
glected, and we use a typical value j3 = 50mm _1 , con- 
sistent with the range of values (3 w 10 — 500mm _1 
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FIG. 4: Diffuse reflectance and the colors Top: color- matching functions x, y, z of the human eye (CIE 1931 XYZ 
standard) as a function of wavelength A (see also Supporting Information). Bottom : calculated diffuse reflectances of CeSF 
and HgS. Shown are the diffuse reflectances of the bulk (Roo) and of a diluted (concentration 20%) thin film of 10/im on a 
white substrate (-Rio M m,20%) computed for /3 =50 mm - . The corresponding color for each case is displayed in the legend). 



reported in the literature for a wide range of indus- 
trial inorganic pigments 35J. For a semi- infinite layer, 
the KM model yields a simple expression for the dif- 
fuse reflectance: ^^(w) = a(oS) — \J a{ui) 2 — 1, with 
a{oS) = 1 + 2cK(u))/f3(u>). As expected, « 1 for 
cKj (3 <C 1 (weak absorption) and R^ ~ for cK//3 1 
(strong absorption). The formula for a sample of finite 
thickness is given in the Supporting Information. 

We have employed the KM model in conjunction with 
our first-principles absorption K(uj) to compute the dif- 
fuse reflectance of CeSF and HgS samples, as well as 
their color (Fig. [3]) - The latter is obtained (e.g. as XYZ 
tristimulus values or xy chromaticities) by taking into 
account the spectral distribution of the light source as 
well as the sensitivity of the human eye to red, green 
and blue light as encoded in the empirical color matching 
functions displayed in the top panel of Fig. @] (for details, 
see Supporting Information). We considered semi- infinite 
samples as well as realistic 10^m layers on a white (i?=l) 
substrate, consisting of either pure (c = 100%) or diluted 
(c = 20%) pigments. The resulting diffuse reflectances 
and colors depicted in Fig. |4] reveal that the pure semi- 
infinite bulk samples of both materials have bright red 
colors, while the thin films with a pigment concentration 
of 20% have a more orange tone, especially for HgS. 

We now introduce a simple performance criterion for 
the usability of such materials as pigments. In order for 
the color to be a bright red, the reflectance for u > A 
should drop sufficiently quickly such that admixtures 
from the green and blue part of the visible spectrum 
are suppressed. This can be insured by requiring that 
i?(A + Mu c )<l/2. Here, hSu c rj 0.2eV is the charac- 
teristic frequency interval over which the human eye dis- 
tinguishes between the primary colors (corresponding to 



5 A s=! 60nm for red light, see top panel of Fig. |4j. Since 
the reflection for lj > A arises from a finite absorption, 
the above requirement can be translated, using the KM 
model, into the criterion: 

K(A + h5u> c ) > /3/(8c) (3) 

This threshold is marked in figure EJb): for pure mate- 
rials (c = 100%) both CeSF and HgS largely satisfy the 
criterion. Indeed, R OQ (X) with /3 = 50mm -1 switches 
from iJoo = to Roo = 1 within only a few nanometers 
(Fig. 0J. Criterion © can be turned, for a given c//3, 
into a minimal thickness L m i n above which the pigment 
has the same apparent color as the bulk, i.e. its color is 
sufficiently stable. This is represented in Fig. [5j Using 
ft = 50 /mm, the minimal thickness of a paint layer con- 
sisting of 20% CeSF is 3.2/im, whereas it is 35/xm for 20% 
HgS. As a result, the reflectivity of 10/im layers with 20% 
diluted HgS violates our quality criterion, which explains 
the notably different (orange) color of those samples, ob- 
served in figure |U 

In conclusion, we presented a framework for the the- 
oretical determination of the color of pigment materials. 
Our methodology combines first-principles calculations 
of the frequency-dependent absorption based on state-of- 
the art many-body theories, with an intermediate length- 
scale modclization of the scattering properties of realistic 
samples. We applied this methodology to classic 'vermil- 
lion' HgS and to the recently discovered rare-earth pig- 
ment CeSF. Our results reveal that the red coloration 
of these two materials has a very different origin at the 
atomic scale. While it is due to large optical matrix ele- 
ments in HgS, the key effect in CeSF is the large density 
of weakly dispersing states (resulting both from the low- 
dimensionality of the crystal structure and the localized 
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FIG. 5: Universal quality curves for thin films of CeSF 
and HgS. The thin films located above the corresponding 
curve in the (/JL m i n ,c//3) coordinates pass the quality crite- 
rion for a good pigment (see the text). The top and right 
axes correspond to a fixed value of j3 of 50/mm, allowing for 
a direct translation of a given concentration into a minimal 
film thickness. 



coupling) , a many-body self-energy is computed and inserted 
into the Green's function of the full solid. Self-consistency 
over the total charge density and the effective atom parame- 
ters is implemented 0, Q. Important technical points of our 
calculational approach are presented in the Supporting Infor- 
mation. 
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character of Ce 4/ states). We also demonstrated that 
the relevant optical transitions in CeSF are S-3p to Ce- 
4/ inter-atomic transitions, in contrast to previous belief. 
Our methodology may lead to rational materials design 
in which first-principles calculations are used in synergy 
with solid-state chemistry in order to create new materi- 
als with specific optical properties. 

METHODS 

Combi ning electronic structure and many-body techniques 
EL Si 0, B l2Q - [22Tj we calculate the electronic structure of rare- 
earth pigments. Then we take these approaches further by 
performing an ab initio calculation of the optical response 
functions and going all the way to a prediction of the actual 
color (expressed e.g. in RGB coordinates [38H39T]) of a sample 
with a given granularity, using an intermediate length-scale 
modelisation [32 1. 

As discussed in the main text, the strong Coulomb repul- 
sion between the 4/ electrons in CeSF is responsible for the 
splitting between the empty Ce-4/ states at the bottom of 
the conduction band and the occupied Ce-4/ states, which 
hybridise with the S-3p states to form the valence band. A 
proper treatment of this strong Coulomb interaction is essen- 
tial in order to correctly reproduce the optical gap and elec- 
tronic structure of this material. Hence, we have employed 
the so-called LDA+DMFT approach [U, [13], which combines 
electronic structure calculations with a many-body treatment 
of local correlations on the Ce 4/ shell in the framework of 
dynamical mean-field theory (DMFT) [2(j. The interaction 
vertex has been computed from first principles using the con- 
straint random phase approximation^, 0]- 

Within LDA+DMFT, the Ce-4/ shell is treated as that 
of an effective atom, self-consistently coupled to an environ- 
ment describing the rest of the solid. From the Hamiltonian 
of this effective atom (which takes into account crystal-field 
effects, intra-atomic Coulomb interactions and the spin-orbit 
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These supporting notes give additional infor- 
mation concerning the calculation of color of 
CeSF and HgS. In particular, we (a) describe 
the technical details of our computational scheme, 
"LDA+DMFT" , the combination of density func- 
tional theory within the local density approxi- 
mation (LDA) with dynamical mean field the- 
ory (DMFT), (b) present the momentum resolved 
electronic structure and optical conductivity of 
our target compounds CeSF and HgS, (c) pro- 
vide further evidence for our interpretation of 
the onset of optical transitions in CeSF being de- 
termined by p-f transitions by discussing the im- 
plications on the electronic structure and optical 
properties that would result from an alternative 
"f-d scenario", (d) summarize the Kubelka-Munk 
theory, which allows us to deduce the color of 
our compounds from the knowledge of the bulk 
optical conductivity, including practical details of 
our intermediate length scale modelisation. We 
conclude with a discussion of the shape of the 
onset of the absorption spectrum for the two 
prototypical compounds that we investigate, the 
quasi-two-dimensional CeSF with its localized f- 
electrons and the three-dimensional, more con- 
ventional HgS with broad and well-hybridising 
bands. 



COMBINED DENSITY FUNCTIONAL AND 
DYNAMICAL MEAN-FIELD SIMULATIONS OF 
CESF 

Our combined density functional and dynamical mean 
field "LDA+DMFT" implementation H Q is based on a 
full-potential linearized augmented plane waves (LAPW) 
technique as implemented in the Wien2k code [3|. Our 
calculations are fully self-consistent, including over the 
charge density jij]. Spin-orbit coupling has been taken 
into account through a second variational treatment. 

The LDA+DMFT calculations for CeSF have been car- 
ried out for the paramagnetic phase at the experimental 
crystal structure [j| ; it is of tetragonal symmetry (space 
group PA/nmm), with the lattice parameters a=3.993A 
and c=6.950A. The correlated / orbitals are constructed 
from projected atomic orbitals which are promoted to 
Wannier function through an orthonormalisation pro- 



cedure. In practice, we have used an energy window 
[—8 : 10.8] eV, which includes bands with mainly / or- 
bital character, following the method implemented for 
d electron materials in Ref. Local correlations be- 
tween the Ce 4/ states have been treated within DMFT. 
The Hubbard and Hund interactions, U = 4.8 eV and 
J = 0.70 eV, resulting from the local Coulomb inter- 
action on the Ce 4/ shell, have been calculated within 
the constrained Random Phase Approximation (cRPA) 
|f|. We have employed the recent implementation of 
cRPA into the LAPW framework of the Wicn2k code 
of Refs. 

The Hubbard-I approximation has been employed as 
a quantum impurity solver within the LDA+DMFT 
scheme. We have used the full four-index Coulomb in- 
teraction matrix in the Hubbard-I impurity solver, thus 
all multiplet effects within the Ce 4/ shell were taken 
into account. Double counting (DC) of Coulomb inter- 
actions between the / electrons was avoided, using the 
fully-localized limit (FLL) expression E CT = U(N— 1/2) — 
J(N/2 — 1/2), where N w 1 is the occupancy of the 
Ce 4/ shell in the DMFT quantum impurity problemQ. 
In the self-consistent LDA+DMFT calculations we em- 
ployed 225 k-points in the irreducible part of the tetrago- 
nal Brillouin zone (BZ) for the BZ-integration. All calcu- 
lations have been carried out at the temperature of 116 K 
with a grid of 1024 Matsubara points, a high-frequency 
expansion up to the tenth order in 1/ (iu> n ) has been em- 
ployed to describe tails of the Green's function and the 
self-energy at Matsubara frequencies above the cutoff. 

The local density approximation leads to a rather 
significant, yet expected, underestimation of the semi- 
conducting gap between S 3p and lanthanide 5d bands 
in the rare-earth fluorosulfides. The LDA+DMFT ap- 
proach employed in this work does not contain correc- 
tions for this LDA error, which is due to long-range in- 
teractions and the missing derivative discontinuity in the 
Kohn-Sham band structure. In addition, CeSF contains 
a partially filled localized 4/ shell, and as pointed out 
in the main text, the relative positions of S 3p, Ce 5d 
and Ce 4/ states define the color of the material. In 
order to evaluate the effect of long-range interactions 
on the relative position of S 3p and Ce 5d bands we 
have performed one-shot calculations within Hcdin's GW 
approximation [Io[ 11 | - u sing the implementation de- 
scribed in detail in Refs. [12|, [l_3| - for CeSF and, for 



sake of comparison, for LaSF. For the latter, the opti- 
cal gap is known to open between the valence S 3p and 
conduction La 5c? bands. The band gap in LaSF calcu- 
lated within the LDA is equal to 1.3 eV compared to the 
experimental optical gap of 2.8 eV @. Within the GW 
approximation, the value of the gap is increased by about 
1.17 eV compared to the LDA one, and, hence, the LDA 
error is compensated by about 80%. As in usual semi- 
conductors, the shape of the highest valence and lowest 
conduction bands, which are crucial for the absorption 
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FIG. 6: The band structure of CeSF in LDA and GW. Band-structure of CeSF obtained within LDA (in red) and GW 
(in green). The Ce 4/ states were treated as valence states in both cases. The bands of Ce 4/ character are shown by dashed 
lines. The topmost S 3p and lowest Ce 5d bands are highlighted by fat lines. One may clearly see that the GW correction 
affects the Ce5d bands much more strongly then S3p ones. 



edge, are barely affected by GW corrections. 

We have performed similar G W simulations fl2 . 13 1 for 
CeSF, and the resulting G W k-rcsolved spectral function 
is compared to the LDA band structure in Fig. [6] As ex- 
pected, the GIF method is not able to describe correctly 
the strongly-correlated Cc 4/ states, which remain metal- 
lic with their bandwidth increased relative to the LDA 
one. The important thing to notice here is that the up- 
ward shift for the bands of Ce 5d character relative to Ep 
is approximately 6 times larger then the corresponding 
(mostly downward) shift of the S 3p bands. Therefore, 
one may conclude that the underestimation of the pd-gap 
in this case is mainly an effect of the missing derivative 
discontinuity in the Kohn-Sham description, and thus a 
too low position of the 5d conduction states. As in the 
case of LaSF, the shape of the topmost valence and lowest 
conduction bands is weakly affected by GW corrections, 
and can be well described as a rigid shift of the corre- 
sponding bands. The value of the gap between those 
states is increased by of about 1.15 eV, similarly to the 
case of LaSF. 

On the basis of our GIF simulations we conclude that, 



in the case of CeSF, corrections to the Kohn-Sham band 
structure of the "uncorrelated" 3p- and 5d-bands and 
thus the relative positions of the CeSF bands relevant 
for its absorption edge - can be well described by a rigid 
shift of the Ce 5d states relative to Ce 4/ and S 3p by 
a value close to the LDA band gap underestimation in 
LaSF. Hence, we have applied a static shift of 1.5 cV to 
the Kohn-Sham eigenvalues for unoccupied states with 
respect to the occupied bands and partially occupied Ce 
4/ states in our LDA+DMFT simulations of CeSF. 



THE LDA+DMFT K-RESOLVED SPECTRAL 
FUNCTION FOR CESF 



The k-resolved spectral function A^oj) encodes the ex- 
citation energies associated with the addition or removal 
of an electron on top of the many-body ground-state $o : 
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FIG. 7: The LDA+DMFT excitation spectrum, k-resolved spectral function (left panel) and the total (black), partial Ce 
4/ (red), Ce 5d (green), S 3p(blue) density of states (right panel) for CeSF 



i4 k (w) = 5I\ $r c k $ o) S{E r - E +uj + fi) 
r 

+ H( $ r 4 1>o) 2 S(E Q -E r +Lo + [i), (4) 



where the sum runs over all excited many-body states 
<£>r with energies Er, ^ is the chemical potential. When 
correlations are absent the many-body states $o and $r 
are just Slater determinants, and the corresponding ex- 
citation energies are one-electron energies of the corre- 
sponding Hartrec-Fock (or Kohn-Sham) problem. 
In the framework of LDA+DMFT and in the basis 
of Kohn-Sham cigenstates the spectral function A^(uj) 
reads: 



A k (w) 



1 



-Im 



Tr (oj + irj + /i - £k - £k(w))~ 



. (5) 

where ek is the diagonal matrix of Kohn-Sham eigen- 
states, £ k (a;) is the DMFT local self-energy upfoldcd 
into the basis of Kohn-Sham eigenstates [l[ (due to this 
upfolding it acquires k-dependence) and 77 is a positive 
infinitesimal. When correlations beyond DFT are not 
included (£ = 0) the k-resolved spectral function has S 
peaks at Kohn-Sham eigenvalues that form one-electron 
Kohn-Sham bands. In the case of metallic systems with 
weak or moderate correlations the Kohn-Sham eigenval- 
ues are shifted and the corresponding S peaks broad- 



ened due to quasi-particle renormalization and life-time 
effects, respectively. In the present case of strongly- 
correlated Ce 4/ states the DMFT self-energy leads to 
the destruction of coherent metallic 4/ states and the for- 
mation of lower and upper Hubbard bands, hence open- 
ing the semiconducting gap. 

The k-resolved spectral function of CeSF along the 
M — T — Z — R path in the Brillouin zone is shown 
in Fig. 1 of the main text. In Fig [7] we show the 
LDA+DMFT spectral function for CeSF for a longer 
(T — X — M — T — Z — R— A — M) path along with 
the corresponding integrated spectral function. 



BAND-STRUCTURE CALCULATION OF HGS 

For cinnabar red, a-HgS, we calculate, as a first step, 
the Kohn-Sham band-structure within the local density 
approximation (LDA), using the Wien2kQ package. The 
crystal structure of a-HgS is hexagonal, with space group 
P3221, 3 formula units per cell, and lattice constants 
a=4.145A and c= 



42 



=9.496A[ 

The resulting band-structure is shown in figure [5] (a) . 
We find an indirect gap of 1.3eV (from valence bands 
near the A-point to the conduction band at T), and an 
optical (=direct) gap of 1.33eV at the T-point. These 
findings are consistent with earlier results using the gen- 
eralized gradient approximation. [IBj Owing to the well- 
established underestimation of insulating gaps within 
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FIG. 8: Crystal and band-structure of a— HgS. (A) crystal structure of a-HgS (mercury atoms in blue, sulfur in yellow, 
graphic produced with Ref [2l[). (B)-(C) Shown are results of (i) DFT-LDA using Wien2k (red, panel (B)) and FPLMTO 
(red, panel (C)), (ii) the G W approximation (green, panel (B)), and (iii) rigidly shifted ("scissors operator") LDA conductions 
bands (blue, panel (B) : Wien2k, panel (C) : FPLMTO). 



density functional theory, these values are much smaller 
than the experimental optical gap of ~ 2.1 eV|l6[. 

In these cases it has become customary (see e.g. the 
recent Ref. [I?}) to employ a scissors operator to tune the 
gap to its experimental value, before proceeding with the 
determination of e.g. optical properties. However, this 
procedure is ad hoc and would deprive us of our claim 
of performing ab initio calculations. Therefore, we went 
a step further and applied the GW approximation [18 
to HgS. GW, a true many-body approach, has proved 
to be a valuable tool for moderately correlated systems, 
and is particularly suited to correct band-gaps in semi- 
conductors 111 11 1 (See also our paragraph on GW ap- 
plied to CeSF in these supporting information). 

For technical reasons, we resorted to a different elec- 
tronic structure technique, the full potential linear muffin 
tin orbital (FPLMTO) method[ig|. Using again the LDA 
potential, the FPLMTO band-structure, shown (in red) 
in figure [8] (b) is in reasonable agreement with our re- 
sults using Wicn2k. Taking these FPLMTO findings as 
starting point, we performed a one-shot G W calculation 



using the code of Ref. [20j . 

The resulting quasi-particle band-structure [as ob- 
tained in the usual way of perturbing the Kohn-Sham 
eigenvalues to linear order in frequency in the GW seli- 
energvfiol]] is shown (as the green dashed lines) in figured] 
(b). The indirect gap comes out to be 2.1eV, in excellent 
agreement with experiment [Tfjj . In the occupied part, 
this GW band-structure closely follows the LDA result 
for energies larger than -2eV. In the figure are also shown 
(in blue) the LDA conduction bands rigidly shifted up- 
wards by a constant value. Evidently, with this scissors 
operation also the conduction bands of the LDA yield a 
favorable approximation to the results over at least 
lcV of the bandwidth of unoccupied states. Therewith 
the "scissored" LDA covers an energy range sufficient to 
account for the onset of optical absorption. 

With this ab initio justification, we turn back to the 
Wien2k data and apply a scissors operator to reproduce 
the GW gap, yielding an optical gap of 2.13eV. The 
resulting shifted conduction bands are shown (again in 
blue) in figure [5] (a). This establishes the band-structure 
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that our calculations of the optical conductivity of HgS 
are based on. 

It is worth mentioning that the onset of optical absorp- 
tion takes place between bands of mainly S3p character 
for the valence band and Hg6s with admixtures of Hg6p 
and S3p for the conduction band. 



OPTICAL CONDUCTIVITIES 



values for the charge gaps for both HgS and CeSF, how- 
ever, we believe this effect to be minor. Of course there 
can be effects of cancellation, e.g. in the case of HgS, be- 
tween the electron- hole interaction induced rcdshift, and 
a potential gap over- estimation within GW. However, 
experience shows that LDA based one-shot GW calcu- 
lations correct well for undcrcstimations of LDA band- 
gaps, and rarely yield a too large gap with respect to 
experiment poj. 



In the main text, we showed the low energy onset of 
absorption in the polarization averaged optical conduc- 
tivity. Owing to the non-cubic crystal structures (CeSF: 
tetragonal, HgS: hexagonal) there is a dependence of the 
optical response on the orientation of the light polariza- 
tion with respect to the crystallographic axes. Figure [9] 
displays our results resolved into the two inequivalent 
orientations (along the crystallographic a and c-axis), 
as well as for an expanded energy range. In both com- 
pounds, the polarization dependence is notable mostly 
at higher energies. In the case of a-HgS we find an 
overall agreement with the dielectric function of Ref. [ItJ 

We further would like to comment on our assumption 
that excitonic effects are negligible in the compounds con- 
sidered, thus verifying our methodology (For a general 
discussion of these effects see Ref. [ll| ) . 

In conventional semiconductors, excitonic effects 
are well-known to severely modify absorption spectra 
through the appearance of sharply defined bound states 
below the conduction bands (see e.g. Ref. [22j for the ex- 
ample of Si02). Such features are absent in the experi- 
mental dielectric function of both HgS[T3|, and CeSFpjj]. 
Besides, the electron-hole attraction in the excitonic 
regime may red-shift the onset of absorption (See e.g. 
Ref. 23J for a discussion on GaAs). Since we obtain good 



F-D INTRA-ATOMIC SCENARIO FOR THE 
ABSORPTION EDGE IN CESF: ELECTRONIC 
STRUCTURE 

Similarly to many other Ce-based semiconductors (for 



example, the sesquioxide Ce203 |24(, or the sesquisul- 
fide Ce2S3 [25-27 1 ) the absorption edge in CeSF has 
been previously ascribed to the Ce intra-atomic 4/ to 5d 
transitions Q. We will call this the "f-d scenario" , in con- 
trast to the true "p-f scenario established by our GW 
calculations for CeSF (see above). As discussed in the 
present manuscript, our LDA+DMFT calculations lead 
thus to a rather different picture of the CeSF electronic 
structure, with the onset of absorption being due to opti- 
cal transitions between the occupied S 3p and empty Ce 
4/ bands, which result in a sharp absorption edge and 
the bright red color of this compound. 

However, we have also simulated the optical proper- 
ties and the coloration for CeSF in the "f-d scenario" by 
using the following procedure: starting from our LDA 
band structure described above, we apply an artificial 
upward shift of 1.5 eV to the partially filled Ce 4/ band. 
This choice is in contradiction to the results of our GW 
simulations, see Sec. of this supporting information, but 
mimicks the situation of other ceria where a gap opens 
between / and d bands. The resulting 11 f-d' electronic 
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FIG. 10: The f-d scenario. The LDA+DMFT band structure of CeSF (left panel) obtained in the f-d scenario. Right panel: 
the corresponding total (black), partial Ce 4/ (red), Ce 5d (green), S 3p(blue) density of states 



structure is depicted in Fig. [TU] 

In this case the occupied Ce 4/ band is located well 
above the top of the S 3p band, hence the absorption 
edge is due to Ce 4/-Ce 5d "intra-atomic" transitions. 
The calculated value of the band gap is 1.52 eV, which 
is smaller than the experimental one. Of course, the gap 
value depends on the choice of the shift. However, the 
sharpness of absorption edge is not sensitive to this pa- 
rameter. One may notice, that, in agreement with previ- 
ous electronic structure calculations (28| . the bottom of 
the conduction band is very dispersive. Therefore, one 
should expect a rather slow growth of the joint density 
of states right above the gap. 

Here, we verify this qualitative analysis and demon- 
strate that the "f-d' scenario fails to account for the es- 
tablished usability of CeSF as a high quality red pigment. 
Shown in hgure [TTJ (a) and (b) are the (polarization av- 
eraged) optical conductivities and absorption coefficients 
for both scenarios. The gaps in both setups differ by sev- 
eral lOOmeV, explaining the different energies at which 
optical absorption sets in. The other major difference in 
the optical spectra is that the intensity in the 11 f-d sce- 
nario" is significantly lower than for the "p-f scenario 
suggested by our GIF calculations (see above). Indeed, 
the absorption coefficient K (lu) shown in figure [TT] (b) 
fails to pass our quality measure for good pigments even 
in the bulk limit (concentration c = 100%) : K(oj) sur- 
passes /3/8 (with /? = 50/mm) at only about 0.2eV above 
the onset. 

The consequences of this large variation in the sharp- 



ness of the absorption onset are displayed in figure [T2J 
which shows the diffuse reflectance and the resulting color 
for both scenarios (the detailed color coordinates are in 
Table 1 : The bulk color (see i?oo) is darker than the 
prominent orange/red of real CeSF that we reproduce in 
the "p-f scenario. This mostly owes to the smaller value 
of the optical gap. As can be seen by looking at the red 
stimulus of the CIE observer, large parts of the red spec- 
trum arc cut off in the p-setup, decreasing the overall 
luminosity. The shallowness of the onset of absorption, 
causes admixtures of 5-10% of green and even blue, which 
slightly increases the brightness, but in a rather uniform 
fashion, i.e. mostly adding components of gray. 



When using "f-d" CeSF as a realistic pigment, i.e. as a 
20% admixture in a layer of 10/xm on a white substrate, 
the color becomes a grayish ocher. This is due to the large 
increase in contributions from all colors, as is evident 
from the diffuse reflectance in figure [12] As expected 
from our analysis (based on our quality measure for the 
absorption coefficient), this instability with regard to the 
deployment in a realistic setup, i.e. the dependence of 
the color on the concentration and the film thickness, 
make "f-d' CeSF unsuitable for practical applications. 
This finding being at variance with the fact that CeSF is 
such a good pigments that it is indeed commercially used, 
thus strengthens our interpretation that in real CeSF the 
onset of absorption is due to p-f transitions, as described 
by our calculations. 



13 




a) [eV] co [eV] 

FIG. 11: The two scenarios of CeSF : impact on the optical conductivity and the absorption coefficient. Com- 
parison of the (true) "p-f the (artificial) "f-cP scenario for (A) the optical conductivity and (B) the absorption coefficient. The 
quality measure /3/8c in the right panel uses again j3 = 50/mm. 
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FIG. 12: The theoretical diffuse reflectance. Comparison of the diffuse reflectance R(X) and the corresponding colors in 
the "p-f and the "/-<?' scenario. 





thickness 


concentration 


X Y Z 


x y 


CeSF 
"p-f scenario 


oo 
oo 
lOfim 


100% 
20% 
20% 


0.60 0.39 0.00 
0.63 0.44 0.01 
0.65 0.47 0.01 


0.60 0.39 
0.58 0.41 
0.57 0.42 


CeSF 
"/-<f scenario 


oo 
oo 
10/im 


100% 
20% 
20% 


0.36 0.26 0.01 
0.55 0.50 0.28 
0.74 0.73 0.44 


0.50 0.37 
0.42 0.37 
0.39 0.38 


HgS 


oo 
oo 
10/im 


100% 
20% 
20% 


0.53 0.33 0.01 
0.60 0.44 0.08 
0.65 0.52 0.08 


0.67 0.33 
0.54 0.39 
0.52 0.42 



TABLE I: The colors of CeSF within the (true) "p-f and the (artificial) "/-(?' scenario, and HgS. Shown are the 
CIE 1931 tristimulus values XYZ and the chromaticities xy (z=l-x-y) of both compounds, for different polarizations, pigment 
concentrations and sample thicknesses on a white substrate. 
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THE COLOR OF MATTER 



The color of a pigment, or any object in general, owes 
to its diffuse, i.e. non-directional, reflectance. This has 
to be contrasted to the specular reflectivity of geometric 
optics as given by Fresnel's formulae. While the latter is 
entirely determined by the absorption properties of the 
system (for a n ap plication to colors and correlated mate- 
rials see Ref. 29|, 



3d]), the diffuse reflectance requires the 
accounting for the scattering of light due to imperfections 
of the specimen, such as impurities, grain boundaries or 
granularities. These are, by definition, effects that are 
beyond a first principles description of extended, perfect 
solids. We therefore resort to a modelisation for objects 
that both absorb and scatter light, where the latter is 
accounted for by an effective scattering rate (3. While 
the latter in principle depends (like the absorption coeffi- 
cient) on the frequency of the light, we shall demonstrate, 
for the case of CeSF, that this dependence is weak and 
non-essential for the fidelity of the calculated color. 



Effective medium theory for light scattering in 
inhomogeneous media — the Kubelka-Munk model 

The model that we employ for our calculations is owing 
to Kubelka and Munk [3lll32| . and is based on considering 
two light fluxes, J + and J_, that account for light prop- 
agating along a coordinate x downwards or upwards in a 
film of thickness L deposited on a substrate of reflectance 
Rs- 

Using the absorption coefficient K = — Im (n + ik) ~ 
and introducing the effective scattering parameter /?, 
the differential fluxes follow (for details see e.g. Refs. [33j, 
H): 



Rn 



a 



vV - 1 



(10) 



dj + {x) 

dx 
dJ_{x) 
dx 



-2KJ + (x)-j3(J + (x) + J-(x)) (6) 
2K J_(x) + (3{J+{x) + J_(x)) (7) 



From this, we obtain for the reflectance R(x) 
— J + (x) / J_(x) at depth x 



-R=^( R+ R- 2 ( 1 + T 



(8) 



integrating x from to L, and thus R from R(0) = 
the substrate reflectance, to R{L), and solving for R(L) 
yields 



R(L) = R^jRs - Roo) - RoojRs - R^)e^ R -- R ^ i{)) 



{Rs - Rc 



{Rs - R^)eP L ( R ™- R ™) 



with a = 1 + 2K/ (3 corresponding to the reflectance of a 
semi-infinite sample (bulk limit). Here, all quantities de- 
pend on the wavelength or frequency of the light. Using a 
pigment as an admixture to a (mostly transparent) host 
material can be modeled, in the simplest possible ansatz, 
by rescaling the absorption with the concentration c, i.e. 
using in the above equations cK instead of K . 



The scattering parameter (3 

The above description on the diffuse reflectance relics 
on the scattering parameter (3 that is not accessible from 
first principle calculations for the perfect solid. In fact, (3 
is a property that depends on the microstructure of the 
material. It arises thus the question of how to determine 
typical values of this parameter in practice, and, in par- 
ticular, whether its frequency dependence is crucial for 
the calculation of the actual color. 

In the case of CeSF there exist experimental measure- 
ments of both, the diffuse reflectance i?[f| and the com- 
plex refractive index n = n + ^fc(28l| . Assuming the valid- 
ity of the KM model, equation (|10|) can thus be inverted 
to extract the scattering coefficient : (3k ai = ~^Tij with 
a expressed as a = \{R + -^). 

Fig. [TBI shows R, K, and the resulting (3km as a func- 
tion of frequency. As it turns out, the frequency depen- 
dence of the scattering coefficient (3 is negligible. Indeed, 
approximating (3 = 50/mm as indicated by the dashed 
line in figure [TBJ the reflectance Rkm within the KM 
model neatly traces the experimental data above the on- 
set of absorption [43l|. This justifies our choice to use the 
approximate (3{cu) = 50mm -1 in all of our calculations. 
This value is moreover consistent with the range of values 
(3 w 10 — 500mm _1 reported in the literature for a wide 
range of industrial inorganic pigments 35 1. 



where 



Opacity and covering power of thin films 

A brilliant color of the bulk, as resulting from R x , 
does not guarantee that a compound can be used as a 
pigment. While in some applications pigments are ad- 
mixtures to the actual material the object is made of, 
pigments are mostly used in lacquer or paint. In other 
words it is not the bulk but the behavior in thin films, 
i.e. the opacity, that defines a good coloring material. 
Thus a good pigment must meet two requirements : 

Imposing the pigment's color. The worst case is a 
"black" substrate : R s = 0, while i?^ ^ 0. Then it 
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FIG. 13: The diffuse reflectance. Experimental diffuse reflectance R of CeSF from Demourges et al.[H], absorption coefficient 
K — w/c Im\/I from the dielectric function e of Goubin et al.[28|], and scattering coefficient Pkm(u), as extracted via the KM 
model. Also shown is the KM diffuse reflectance Rkm that results from the experimental absorption coefficient K and an 
energy-independent scattering coefficient /3 = 50/mm. The magnitude of the latter is indicated by the dashed line. 



follows from equation 



R(L) = R c 



e /3L(R^ 1 —Rr, 



R 2 



< Ro 



(11) 



As a matter of fact the finite thickness of the film results, 
for the current case, in a more sharply defined absorp- 
tion : While R(L) = 1 for R x = 1, R(L) < Roo for 
Roo < 1. 

Hiding the color of the substrate. Now the worst case 
is a "white" substrate : Rs = 1, while R^ <C 1. Then 



R(L) = 



l + i?ooe 



Rc 



(12) 



In that case R(L) > i?oo, which means that a finite film 
will have a broader tail in the reflectance and admixtures 
of wavelengths away from the absorption edge arc larger 
than for the bulk material. Hence the quality measure 
of a pigment that we introduce in the main manuscript 
considers a white substrate. 



From the diffuse reflectance to the perceived color 

While many spectroscopic devices have been devised to 
analyse matter, the spectrometer most commonly used is 
the human eye : it sensors the response to incident light 
in a range from roughly A = 400 to 700 nm, and the 
information gathered by the rods and cones of the retina 
are post-processed by our brain into a brightness and 
color sensation. 

The calibration, or color matching functions, of the 
three types of cones, responsible for the color perception, 
can be empirically measured (3fj| . and used to setup 
a three-dimensional color space. Many variants of color 



schemes exist, and they are basically related by a unitary 
transformations of their coordinates. In the pre sent work, 
we employ the CIE 1931 XYZ color space [31, Hi- The 
corresponding matching functions, probing the response 
to red, green, and blue light, are denoted x(A), y(A) and 
z(A), and are displayed in the top panel of Fig. 4 of the 
manuscript or, here, figure [T2l 

Knowing the reflectivity R(X) of a material or object 
(as e.g. in the bottom panel of Fig. 4 of the main text 
or, here, figure IT2"j) . and the spectral distribution iS'(A) of 
the light source, i.e. accounting for the type of light that 
the considered object is illuminated with, the XYZ color 
coordinates, or tristimulus values, are calculated as (see 
e.g. Ref. Ei|) 



X,Y,Z) = k / d\R(X)S(X) 5(A),y(A),z(A) 13) 



where k is a normalization constant chosen such that 
k J dXS(X)y(X) = 1. As light source we use in our cal- 
culations the CIE standard illuminant D65 [4l|, which 
is appropriate for usual day light. While the tristimulus 
coordinates include information on both color and bright- 
ness, in color pictograms often only the color information 
given by the normalized chromaticity values 
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(14) 



which arc shown in Table 1 for the two scenarios of CeSF, 
as well as for HgS. 
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INTERBAND ABSORPTION IN TWO-BAND 
MODEL 



Let us now introduce the band dispersions corresponding 
to the "spherical" (Fig. Hih ) 



In order to understand how the topology of the bottom 
of conduction/top of valence bands in CcSF influences 
its optical properties we consider a simplified model, in 
which we include a single valence and a single conduction 
band with analytical dispersion laws. This model is ade- 
quate to describe the optical conductivity in the vicinity 
of the absorption edge, where the relevant interband tran- 
sitions are between single valence and conduction bands, 
which have a quasi-two-dimensional shape in the case of 
CeSF (see Fig. [Hp)- Thus, to model the absorption edge 
in CeSF we consider the case of "cylindrical" valence and 
conduction bands (Fig.fMk). Three-dimensional ("spher- 
ical") bands (Fig.Q3b) provide a more adequate model of 
the actual states involved in optical transitions at the ab- 
sorption edge of HgS (especially for the conduction band, 
see Fig. EH) . 

In calculations of the interband absorption one starts 
from the real part of the optical conductivity (Eq. 1 in 
the main text), which in the case of zero temperature can 
be written as: 

Reo-M = V / — tr[v(k)A(k,Lj)v(k)A(k,w+v)]. 

V J v 

(15) 

The imaginary part of the dielectric function (in CGS 
units) and absorption K are then evaluated from Reer(cj) 
as 

Ime^( M ) = — R C g c 9-» = Rea ^ (1 6 ) 



K{u>) 



uik(uj) wIme C9S (w) Rea(w) 
c 2nc 2nceo 



(17) 



where h = n + ik is the refractive index, c is the speed 
of light, eo is the vacuum permittivity. The absorption 
K(uj) has the dimension of inverse length. 

Neglecting the k-dependence for the velocities, 
W12 (k)«2i (k) = — v 2 one may carry out the energy in- 
tegration in (|15[) and obtain the optical conductivity ex- 
pressed through the "joint density of states" : 

J dud [hw - E v (k)] 6 [thj - E c {k) + hv] 

= ls[E v (k)-E c (k) + hu], (18) 

where E v (k) and E v (k) are the valence and conduction 
bands dispersions, respectively. Then (fl~5j) reads as fol- 
lows: 



Recr(^) = / *t_* s [Ev{k) _ Ec{k) + H 



E v (k) = - 



h 2 k 2 

2m v 



E c (k) = t^- 4 



2m r 



and "cylindrical" (Fig. [Llk 
E v (k) = - 



h 2 {k 2 x + k 2 y ) 
2m v 



E c (k) 



h 2 (k 2 



2m r 



(20) 



(21) 



(22) 



(23) 



cases. Here, m v r c \ are the valence (conduction) band's 
effective masses and A is the band gap value. 

By substituting the d 3 k integration in (|19p with the in- 
tegration over the constant energy surface / d 3 kS[E(k) — 
E ]... — s- J v fc [gfif)-g ] ••• one obtains the real part of the 
optical conductivity: 



Rea sph (Su 



2V2e 2 v 2 

1 

(A + h5cu)irh 2 



.3/2 



Rea cyl (5. 



2e 2 v 



2„,2 



(A + hSuj)ah 



(HSlo) 



1/2 



(24) 



(25) 



for the spherical and cylindrical dispersions, respectively. 
Here, m r is the reduced effective mass m r = m c m v /(m c + 
m v ). hSu is the photon energy relative to the absorption 
edge, Hv = A + Mw. Of course, Rc<j(8lu) = within 
the gap (Su> < 0), hence one may immediately see that 
Rea cvl (5iv) exhibit a jump at the absorption edge, while 
Rea- sph (Su>) is continuous everywhere. 

The corresponding absorptions (fT7| for Slj > (for 
8bj < the absorption is again zero) read 



K sph (5u>) 



V2e 2 v 2 



7r(A + hSttj)hh 2 ceo 



3/2 



(HSu> 



1/2 



K cy i(6uj) 



e 2 v 2 



(A + h8uj)ahnceQ 



m r . 



(26) 



(27) 



(19) 



To model the CeSF case one may assume a flat con- 
duction band (the dispersion of the Ce 4/ states at the 
bottom of conduction band is negligible, see Fig. H4b). 
hence m r = m v . For the effective mass we used the 
value m r = m v equal to 0.6too, which was extracted 
from the curvature of the topmost S 3p band at T in the 
CcSF LDA+DMFT momentum-resolved spectral func- 
tion, Fig. [7J Other parameters entering in Eqs. (|26I27[) 
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FIG. 14: Upper panel: a. Two model bands "cylindrical" along the Y — Z path (along other two paths Y — X and Y — Y 
they have the same shape as shown in c).; b. Zoom at the CeSF bands along the Y — Z path.; c. Two model bands in the 
spherical case; d. Zoom at the HgS bands along the M — Y — A path. Lower panel: Absorption K(Soj) for the two-band models 
with spherical and cylindrical band dispersions. The cyan and pink horizontal lines indicate the threshold value of K for the 
semi-infinite bulk and thin film cases, see the text 



were chosen as follows: A=2.1 eV, h rj 2 was ob- 
tained from the full complex optical conductivity of 
CeSF, the average Fermi velocity can be estimated to 
be v = V p f ■ 5ce~s ~ 0.3 eVA, where 5c e -s ~ 3A is the 
distance between next-nearest Ce and S sites in CeSF, 
and Vpf is a hybridization matrix element between the 
involved orbitals. 

The resulting absorptions for both the spherical and 
cylindrical cases are displayed in Fig. [H] with the same 
set of model parameters employed for both cases. One 
may see that the stepwise behavior of the curve in the 
"cylindrical" case results in a much higher absorption 
just above the edge as compared to the spherical case. 
This leads to an important difference in the color prop- 
erties in the two cases. To quantify those properties we 
again employ the KM model with the scattering param- 
eter P = 50/mm and compute the characteristic absorp- 
tion thresholds for the two cases considered in the article: 
the semi-infinite bulk of pure pigment and a thin layer of 
10 /Ltm with the pigment concentration c =0.20. In the 
"spherical" case the characteristic absorption threshold 



<5wi/2, for which the diffuse reflectance R = 1/2, is equal 
to ~7 meV for the semi-infinite bulk. Hence, the absorp- 
tion edge is sufficiently sharp in this case. However, in 
the case of a thin film we obtain Suji/2=0.52 eV, which is 
substantially larger than the characteristic frequency in- 
terval 8lo c «0.2eV (see the main text). Therefore, the 
"model material" with the spherical bands (Fig. IT4"b ) 
would hardly make a very good pigment. In contrast, in 
the model system with the two-dimensional band topol- 
ogy the diffuse reflectance will undergo a discontinuous 
transition from R = 1 to to R w exactly at the edge. 
Hence, a quasi two-dimensional topology of the conduc- 
tion/valence bands naturally leads to a sharp onset of ab- 
sorption and drop in R as predicted by the LDA+DMFT 
simulations for CeSF (see Fig. 4 in the main manuscript). 
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